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Abstract. Bound states of the Hamiltonian describing a quantum particle living on 
three dimensional straight strip of width d are investigated. We impose the Neumann 
boundary condition on the two concentric windows of the radii a and b located on 
the opposite walls and the Dirichlet boundary condition on the remaining part of the 
boundary of the strip. We prove that such a system exhibits discrete eigenvalues 
below the essential spectrum for any a, b > 0. When a and b tend to the infinity, the 
asymptotic of the eigenvalue is derived. A comparative analysis with the one-window 
case reveals that due to the additional possibility of the regulating energy spectrum 
the anticrossing structure builds up as a function of the inner radius with its sharpness 
increasing for the larger outer radius. Mathematical and physical interpretation of 
the obtained results is presented; namely, it is derived that the anticrossings are 
accompanied by the drastic changes of the wave function localization. Parallels are 
drawn to the other structures exhibiting similar phenomena; in particular, it is proved 
that, contrary to the two-dimensional geometry, at the critical Neumann radii true 
bound states exist. 
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1. Introduction 

The task of finding eigenenergies E n and corresponding eigenfunctions / n (r), n— 1,2,... 
of the Laplacian in the two- (2D) and three-dimensional (3D) domain Q with mixed 
Dirichlet 

/»(r)l«i D = (1) 

and Neumann 

nV/ n (r)| 9njv = (2) 

boundary conditions on its confining surface (for 3D) or line (for 2D) d£l = d^o^d^N (n 
is a unit normal vector to d£l) [1-12] is commonly referred to as Zaremba problem [13]. 
Apart from the purely mathematical interest, an analysis of such solutions is of a large 
practical significance as they describe miscellaneous physical systems. For example, the 
temperature T of the solid ball floating in the icewater obeys the Neumann condition 
on the part of the boundary which is in the air while the underwater section of the body 
imposes on T the Dirichlet demand [3]. Mixed boundary conditions were applied for 
the study of the spectral properties of the quantized barrier billiards [5] and of the ray 
splitting in a variety of physical situations [1]. The problem of the Neumann disc in 
the Dirichlet plane emerges naturally in electrostatics [14]. In the limit of the vanishing 
Dirichlet part of the border the reciprocal of the first eigenvalue describes the mean 
first passage time of Brownian motion to 8VLd [12]. In cellular biology, the study of 
the diffusive motion of ions or molecules in neurobiological microstructures essentially 
employs the combination of these two types of the boundary coniditons on the different 
parts of the confinement [15]. 

One class of Zaremba geometries that recently received a lot of attention from 
mathematicians and physicists are 2D and 3D straight and bent quantum wave guides 
[16-41]. In particular, the conditions for the existence of the bound states and resonances 
in such classically unbound system were considered for the miscellaneous permutations 
of the Dirichlet and Neumann domains [21,26,34]. Bound states lying below the 
essential spectrum of the corresponding straight part were predicted to exist for the 
curved 2D channel if its inner and outer interfaces support the Dirichlet and Neumann 
requirements, respectively, and not for the opposite configuration [23,25,28]. This was 
an extension of the previous theoretical studies of the existence of the bound states 
for the pure Dirichlet bent wave guide [42,43] that were confirmed experimentally 
[44]. Magnetic field influence on the Dirichlet-Neumann structures was analyzed too 
[29,34,36]. Also, for the 2D straight Dirichlet wave guide the existence of the bound 
state below the essential spectrum was predicted when the Neumann window is placed 
on its confining surface [16,18]. From practical point of view, such configuration can 
be realized in the form of the two window-coupled semiconductor channels of equal 
widths [16, 19] whose experimental creation and study has been made possible [45] 
due to the advances of the modern growth nanotechnologies. The number of the 
bound states increases with the window length L and their energies are monotonically 
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decreasing functions of L [24]. Recently, this result was extended to the case of the 
3D spatial Dirichlet duct with circular Neumann disc [39] for which a proof of the 
bound state existence was confirmed, the number of discrete eigenvalues as a function 
of the disc radius a was evaluated and their asymptotics for the large a was given. As 
mentioned above, such Zaremba configuration is indispensable for the investigation of 
the electrostatic phenomena [14]. Similar to the 2D case, it can be also considered as the 
equal widths limit of the two 3D coupled Dirichlet ducts of, in general, different widths 
with the window in their common boundary [19, 32]. Another motivation stems from the 
phenomenological Ginzburg-Landau theory of superconductivity [46] which states that 
the boundary condition for the order parameter \I>(r) of the superconducting electrons 
reads 



n being an inward unit vector normal to the confining interface dQ. Extrapolation 
length (or de Gennes distance) A is equal to infinity for the superconductor/dielectric 
boundary and tends to zero for the contact with ferromagnets [46]. Thus, placing atop 
the superconductor the materials with these two limiting extrapolation lnegths, one 
inevitably needs to deal with mixed boundary conditions. 

In the present research, we discuss the case of the two concentric Neumann windows 
with, in general, different radii a and b on the opposite walls of the 3D straight wave 
guide. A comparative analysis with the one Neumann disk reveals that the two- 
window geometry offers more possibilities for the varying spectral and localization 
properties by the additional channel of changing the second disc radius. Namely, its 
variation leads to the formation of the anticrossing structure of the energy spectrum as 
a function of smaller radius b when the almost flat parts of the energy 6-dependence are 
alternated by their sharp drops down to the lower level pushing it to occupy the next 
lying below neighbouring state, etc. Transitions get sharper and the gaps between the 
anticrossing states decrease with increasing the outer radius a. These avoided crossings 
are accompanied by the drastic changes of the localization of the wave function. 

The rest of the paper is organized as follows. In section 2, we define the model and 
recall some known results. In section 3, we present the main result of this investigation 
followed by a discussion. Section 4 is devoted for numerical experiments, and concluding 
remarks are collected in section 5. 

2. Model and Formulation 

The system we are going to study is given in figure 1. We consider a Schrodinger particle 
with mass m p whose motion is confined to a pair of parallel planes separated by the 
width d. For simplicity, we assume that they are placed at z = and z = d. We shall 
denote this configuration space by Q 



nVf 




(3) 



Q = R 2 x [0,d]. 



Di 



4 



Dirichlet boundary condition 




Figure 1. Dirichlet wave guide with two concentric Neumann disc windows on the 
opposite walls with (in general) different radii a and b. 

Let 70(a) be a disc of radius a with its center at (0,0,0) and 7<j(6) be a disc of 
radius b centered at (0, 0, d), where 

1 .(z) = {(x,y,-)eR 3 ; x 2 + y 2 <z 2 }. (4) 

Without loss of generality we assume that < b < a. We set T = dfl\(j (a) U 7<j(&))- 
We consider Dirichlet boundary condition on T and Neumann boundary condition in 
7o(a) and 7<j(6) what means that <9f2£< and cft|/v from (1) and (2) take the following 
form: dVl D = T, d£l N = 70(a) U 7d(6). 

S.i. TTie Hamiltonian 

Let us define the self-adjoint operator on L 2 (f2) corresponding to the particle 
Hamiltonian H. This is done by means of the quadratic forms. Namely, let go be 
the quadratic form 

qo(f\g)= I W'Vffdr, with domain Q(q ) = {f G H\Q); f\dtt D = 0}, (5) 
Jn 

where if 1 (fi) = {/ G L 2 (fi)|V/ G L 2 (f2)} is the standard Sobolev space and we denote 
by f\dQ,£> the trace of the function / on dQo- It follows that go is a densely defined, 
symmetric, positive and closed quadratic form. We denote the unique self-adjoint 
operator associated with g by H and its domain by D(Q). It is the Hamiltonian 
describing our system. From [47] (page 276), we infer that the domain D(Q) of H is 

D(Q) = {fe H\n); -Af G L 2 (Q),f\dQ D = 0,^\dQ N = o} 
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and 

Hf = -Af, V/GD(fi), 
where we have set h 2 /(2m p ) = 1. 

2.2. Some known facts 

Let us start this subsection by recalling that in the particular case when a = b = 0, 
we get H°, the Dirichlet Laplacian, and at a = b = +oo we get H°°, the Neumann 
Laplacian. Since 

H = (-A R2 ) <g> I © I <g> (-A M ) , on L 2 (M 2 ) <g> L 2 ([0, d]), 

(see [47]) we get that the spectrum of H° is (^) 2 ,+oo and the spectrum of H°° is 
[0, +oo[. Consequently, we have 

Q)\+oo C a (if) C [0,+oo[. 

This leads to a natural choice of the unit of energy as (ir/d) 2 used below. In addition, 
if not stated otherwise, we will measure all distances in units of the wave guide width 
d, and all momenta, in units of 1/d. 

Using the property that the essential spectra are preserved under compact 
perturbation, we deduce that the essential spectrum of H for any finite a and b is 

a ess (jl^j = [l,+oo[. 

An immediate consequence is that the discrete spectrum, if it exists, lies in [0,1]. 

2.3. Preliminary: Cylindrical coordinates 

As the system has a cylindrical symmetry, it is natural to consider the cylindrical 
coordinates system r = (r,9,z). Indeed, we have that 

L 2 {Q, dxdydz) = L 2 (]0, +oo[x [0, 2tt[x [0, l],rdrd6dz). 

We denote by (, ) r the scalar product in L 2 (Q, dxdydz) = L 2 (}0, +oo[x [0, 2ir[x [0, 1], rdrdOdz) 
given by 



(f,g)r = / fgrdrdOdz. 

«/10,+oo[x[0,27r[x[0,ll 



]0,+oo[x[0,2tt[x[0,1] 

The Laplacian operator is written as 

_ 2 _ ld_ f d_\ ]_cP_ 
rAz ~ r ~ rdr\dr) + r* dQ 2 + dz* ' { } 

Therefore, the eigenvalue equation is given by 

-A r , diZ f(r,9,z) = n 2 Ef(r,9,z). (7) 
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Since the operator is positive, we set ti 2 E = k 2 . Equation (7) is solved by separating 
the variables and considering the function / as a product 

f(r,9,z) = R(r)-&(9)-Z(z). (8) 

Plugging the last expression into (7), one first separates Z by putting all the z 
dependence into one term so that ^ can only be constant taken, for convenience, as — s 2 
with its value determined by the boundary conditions at z = and z — 1. Second, we 
separate the term ^- which has all the 9 dependence. Using the fact that the problem 
has an axial symmetry and the solution has to be 2n periodic and single-valued in 9, 
we obtain that ^- should be a constant — m 2 for m E Z. Finally, we get the following 
equation for R 

R"(r) + l#(r) + (V - s 2 - ^pj R(r) = 0. (9) 

We notice that equation (9) is the Bessel equation and its solutions could be expressed 
in terms of the Bessel functions [48,49]. More explicit solutions could be given by 
considering boundary conditions; for example, for the Dirichlet requirements at the 
both walls one has sf D = (j + l)n while for the NN situation sf N = jir, and for the 
mixed case s^ D = (j + l/2)7r, j E N (see Section 4 for more discussion). 

3. Analytical results 

Here we prove existence conditions and provide evaluations derived from the analytical 
consideration. 

Theorem 3.1 The operator H has at least one isolated eigenvalue in [0, 1] for any 
nonzero a and b. 

Moreover, forb big enough and X(a,b) being an eigenvalue of H less then \, there 
exist positive constants C a and C& such that 

AM), (%,%). (10) 

Remark 3.2 The first claim of Theorem 3.1 is valid for more general shape of bounded 
surface S with Neumann boundary condition, not necessarily a disc; it suffices that the 
surface contains a disc of radius a > 0. 

For the more general shape S, using discs of radii a and a' such that 

7o (a) C S C 70 (a'), (11) 

and a comparison argument, one gets the localization of the discrete spectrum (see also 
Refs. [19,32]). 

Proof. For the proof of the the first claim one may mimic the argument given in 
[39] and adjusted for the case of the two windows; however, much simpler and elegant 
way is to use the fact that the Neumann window is a negative perturbation [19]. Thus, 
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if the one Neumann window creates the bound state [39], the insertion of the second 
one just pushes it lower. ■ 

The proof of the second claim is based on the bracketing argument. Let us split 
L 2 (n,rdrd6dz) as follows: L 2 (n,rdrd6dz) = L 2 (tt- b ,rdrd6dz)®L 2 (tt+ b ,rdrd6dz), with 

Q~ b = {(r,e,z) e [0,b+(a-b)(l-z)} x [0, 2tt[x [0, 1]} , 

^ = 

Therefore, 

Here, we index by D and N depending on the boundary conditions considered on the 
surface r = b + (a — b)(l — z). The min-max principle leads to 

a ess (h) = a ess (if+f ) = a ess (h^) = [l,+oo[. 

Hence, if H~' b D exhibits a discrete spectrum below 1, then H does as well. We mention 
that this is not a necessary condition. If we denote by Xj (^Ha, b D ^j > (-^0^) anc ^ 

Xj (jl^ t ne J _ th eigenvalue of H^' b D , H~' b N and H, respectively, then, again, the minimax 
principle yields the following 

A, (H~ b N ) < A,- (H) < Aj (^ 6 D ) (12) 

and for j > 2 

Vi < Aj (h) < Aj (jf"^) . (13) 

As the computation of the eigenvalue of a frustum is not an easy task, let us remark 
that 

Aj (H-f) < X, (H7f) < X, [H b f) (14) 

(the same chain of inequalities is true for the corresponding Neumann Hamiltonians 
H~' N , too). Then, from equations (13) and (14) it follows: 

Vi {Kf) < Aj (h) < Aj ( V) • (15) 
The Hamiltonian -f^ has a sequence of eigenvalues [48, 49] given by 

2 



, \2 , I x \m\n 
Xlmnii) = [In) + 



where x\ m \ n is the nth positive zero of Bessel function of the order \m\ [48, 49] and index 
j amalgamates the three quantum numbers: transverse /, radial n and azimuthal m 
ones, j = (l,m,n). The condition 

kmn < 1 (16) 
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yields that I = 0, so we get A 0m „(&) = (x\ m \ n /b) 2 and A 0m ' n '( a ) = {x\m>\n'/af with 
primed n and m referring to the index j — 1 in (15). Accordingly, for any eigenvalue 
\(a,b) of the Hamiltonian H, there exist m,n,m',n' € N such that 

^<A(a,6)<^. (17) 

When b is big enough, we get the result. ■ 

The above derivation shows that the coefficients C a and C& in equation (10) are, 
actually, the squares of the zeros of the Bessel functions. 

Remark 3.3 As it is shown in the next section, for the large enough diameters 2b many 
bound states exist which are characterized by the azimuthal m and radial n quantum 
numbers. For each of them the estimate (10) holds; however, the magnitude of the 
smaller radius above which this estimate becomes true is different for the different states 
(n,m). Equation (17) and the properties of the zeros of the Bessel function [48,49] 
manifest that estimate (10) for the higher lying states is achieved at the larger b. 

Next statement describes the case when the new bound state just emerges from the 
continuum what means that its energy is equal to the Dirichlet fundamental propagation 
threshold. We assume that the inner radius is either fixed in the interval < b < a or 
is equal to its outer counterpart a. 

Proposition 3.4 When the radius a is equal to a critical value ai at which a new bound 
state emerges from the continuum, equation (7) with E = 1 has a solution fj-°\r,9,z), 
unique to a nonzero multiplicative constant which at infinity behaves like 



p iraO 

/, V >,M) = r 



; ( )/ - . e 



+ p l — 7 ^sm2nz + O 



- - v , ..... , . /• > x. (18) 

with some constants fy. Here, the two quantum numbers n and m are compacted into 
the single index I: I = (n, m). 

A proof of this statement will be given in the next section. 

Note that the derived asymptotics from (18) drastically differs from its quasi-one- 
dimensional couterpart [16, 17,24,26]; namely, if, in the latter case, the state at E — 1 
is always a resonance, for the 3D geometry considered in this paper this is true for the 
azimuthal quantum number m = only, and for \m\ > 2 one has true bound states with 
square integrable functions f^ (r). We point out that the function of the state with 
\m\ — 1 vanishes at infinity, however, it is not square integrable over the whole r — 
plane; accordingly, one can consider it as a resonance too. In other words, a behaviour 
of the states lying on the border between the essential and discrete spectra crucially 
depends on the dimensionality of the problem. 



4. Numerical computations: mode-matching method 



This section is devoted to some numerical computations. We analyse a dependence of 
the eigenenergies E and coresponding eigenfunctions / on the radii a and b in the whole 
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range of their variation. We also compare the obtained results with the one-window 
case [39]. 

4-1. Two identical Neumann windows 

Consider first the case of the equal radii on the upper and lower walls of the duct. The 
eigenvalue equation (7) uses the Laplacian provided in (6). A general solution is written 
in equation (8). Since for our problem a cylindrical symmetry is conserved, one has 
0(0) = -j=e im6 . A total solution can be written: 



(19a) 



/2w 

for one Neumann window: 



f<(r,6,z) = ;^e^f>] m| J H - (j + l/2) 2 r) x? D (z), r < a 



3=0 

00 



Mr, 6, z) = -±=<T» g A™K M (vyj(j + l) 2 - Er^j X f D (z), r > a; (19b) 

for the two identical concentric Neumann windows of the radius a on the opposite 
walls: 

1 00 

/<(r,M) = 4=^^SHJ H (vrv^fr) X f N (z), r < a (20) 



2tt 

3=0 

and for r > a it coincides with (196). Here, J m (r) is Bessel function of the first 
kind with K m (r) being its modified counterpart of the second kind [48,49], coefficients 
and present a relative contribution of the corresponding component into 
the total wave function, and the orthonormalized functions xf D (. z )i xf°( z ) an d 
xf N (z) describing a transverse motion in the Dirichlet-Dirichlet, Dirichlet-Neumann 
and Neumann-Neumann strips, respectively, read: 

X f D (z) =V2sm(j + l)7rz (21a) 
X f D (z) = V2 cos ( j + l/2)nz (216) 



Their corresponding eigenvalues are, respectively, Ef D = (j + l) 2 , Ef D = (j + 1/2) 2 
and Ef N = j 2 . Note that for j ^ the Bessel functions J H E - (j + l/2) 2 r^ 

in (19a) and J\ m \ E — j 2 r j in (20) transform into their modified counterparts 

J H ^(j + l/2) 2 -Er^ and J H (iry/f - Er^j [48,49], respectively. 

Physically, the emergence of the bound state for at least one nonzero radius is due 
to the fact that the introduction of one or two Neumann windows mixes longitudinal 
and transverse motions in the wave guide what results in the splitting off of the unbound 
level from the lowest subband down below the fundamental propagation threshold and 
its corresponding transformation into the localized state with square-integrable wave 
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function vanishing at infinity. Mathematically, this mixing is reflected in the fact 
that the summation index j in equations (19) and (20) enters each component of both 
transverse Z(z) as well as radial R(r) parts of the wave function binding them together. 
Note that this binding is different inside and outside the disc. In a sense, a Neumann 
perturbation of the Dirichlet duct presents a shelter where the Schrodinger particle 
can dwell with its momentum smaller than the cutoff momentum of the fundamental 
propagation threshold. Similar to the bend in curved Dirichlet wave guide [50], one can 
draw parallels with the creation by the obstacle of an attractive quantum well that binds 
the particle; namely, the wave can still be propagating inside the circle but vanishing 
outside: as equations (19a) and (20) demonstrate, there is one propagating along the 
radial direction component [the term with j = when the difference E — (j + 1/2) 2 (for 
one window) or E — j 2 (for the two discs) is positive] while outside the discs the wave 
function possesses exponentially vanishing contributions only. 

Matching function / and its radial derivative at r = a, one gets: 



j'=0 



(22) 



with the following matrix elements of the infinite square matrix Q: 



Qjji — 



J\m\ ( *\/E - (f + l/2) z a 



(3 + 1) 



E 



K \M\«VV + 1 ) ~ Ea 



J H (ny/E-(j' + l/2) 2 a?) 
K\ m \ (ny/{j + l) 2 -Ea^ 



p(i) 



(23) 



where the prime denotes a derivative of the Bessel function with respect to its argument 
and the superscript index i — 1(2) corresponds to one(two) Neumann window(s). The 
infinite square matrices and P*- 2 * 1 present the quantitative measure of the coupling 
between different modes of the DD and ND [for P^] or DD and NN [for p( 2 )] channels: 



P (i) 



p(2) 
r 33' 



f D (z) X f, D (z)dz 

NN i 



3 + 1 



7T (j + l) 2 - (f + 1/2) = 



(24) 



{z) X » N {z)dz 



3+1 



(3+1) 2 



f = o. 



0, jV0,j' = j + 1 (25) 

Nontrivial solution of system (22) exists only when the determinant of the matrix 
Q is an identical zero: 

det||Q|| = 0. (26) 
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This is an equation for determining the bound-state energies E for the different radii a. 

From (22) the coefficients can be determined and, next, from the matching 
conditions, the values of are calculated too. Normalization condition 

/ dO dz drr\f(r,9,z)\ 2 = 1 (27) 
Jo Jo Jo 

allows one to fully construct the function /. Note that the radial integrals in (27) with 
the function / given by equations (19) and (20) are calculated analytically [51,52] in 
terms of the Bessel functions themself. 

Fig. 2 depicts bound-state energies as a function of the normalized radius a. Both 
configurations show a lot of similarities. As discussed in [39] for one and earlier in this 
paper for the two identical Neumann windows, at least one bound state exists for the 
arbitrarily small radius a: in both cases this is the state with the azimuthal quantum 
number m — 0. For quite small radii its energy lies very close to the fundamental 
propagation threshold of the uniform Dirichlet wave guide Eq D — 1. It was conjectured 
[19] that for the one window the energy of this single state at the small radius depends 
on a as 

E(a) = 1 -exp (-^) , a^O. (28) 

Our results show that the constant c is in the range 0.44 < c < 0.45. Large numerical 
errors for the energies close to the threshold precluded more accurate evaluation. 

Further growth of a forces the bound-state energy to decrease, and at large radii it 
tends from above to the fundamental propagation threshold of the Neumann-Dirichlet 
Eq D = 1/4 (for one window, ref. [39]) or the uniform Neumann Eq N = (for two 
windows, see above) channels. Second bound state possesses quantum number \m\ — 1 
and emerges from the continuum at a ~ 1.027 for one window and a ~ 0.866 for the two 
Neumann discs. This difference is not surprising since the two Neumann regions present 
a stronger perturbation to the particle motion in the uniform Dirichlet wave guide as 
compared to the one window; accordingly, in the former case the second bound state is 
formed at the smaller radius. This is true for all other states (but the first one) as well; 
for example, the third bound state with \m\ — 2 emerges for the one disc at a ~ 1.549 
and for the two Neumann windows - at a ~ 1.319. As our results show, in both cases 
the next bound level is another state with azimuthal quantum number m = preceding 
in this way its counterpart with \m\ = 3. Similar to the ground state, energies of the 
higher lying levels decrease when the radius increases with the difference between them 
diminishing what means the increase of the corresponding density of states (number of 
states per unit energy) M{E) until at a = oo one arrives at the continuous spectrum of 
the ND (one window) or NN (two windows) channel with the diverging density M{E). 

Figures 3 and 4 depict the lowest bound-state wave function f(r,z) for one and 
two windows, respectively, at several Neumann radii a. Since the lowest level possesses 
azimuthal quantum number m — 0, its wave function is ^-independent. It is seen that 
for the small radii, say, a = 0.3 in figures 3 and 4, the wave function is almost flat 




Figure 2. Bound-state energies E as a function of the radius a for (a) one Neumann 
window and (b) two concentric Neumann discs on the opposite sides of the wave guide 
where solid lines are for the states with the azimuthal quantum number m = 0, dotted 
lines - for the states with \m\ — 1, dashed lines - for \m\ = 2, and dash-dotted curves 
- for 1 77i I = 3. Note different energy scales in panels (a) and (b). 
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Figure 3. Ground bound-state wave function / in terms of r and z for one Neumann 
window and several values of the radius a. 

spreading far away from the Neumann area. This can be easily understood by means of 
the mentioned above analogy with the quantum well [53]; namely, for the small radius 
a the effective positive potential created by the disc, only barely binds the particle with 
its wave function relatively slowly vanishing at infinity. Increasing Neumann radius 
strengthens the attractiveness of the well; accordingly, as the panels for a = 0.4 (figure 
4) and a = 0.5 (figure 3) exhibit, the function / is stronger localized inside the disc 
with the magnitude of its maximum growing. In both configurations the maximum is 
achieved at the centre of the disc, r = 0. As our calculations show, for one window 
this extremum reaches it maximal value of f max ~ 1-28 at a ~ 0.74 and for the two 
discs - at a ~ 0.56 with f max ~ 1.12 (see corresponding panels in figures 3 and 4). 
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Subsequent increase of the Neumann area makes the potential well wider and wider 
with the wave function / being almost uniformly distributed inside it. Its radial slope 
diminishes, as a transition from the situation with a = 1 to the one with a = 3 and, 
next, to a = 5 in figures 3 and 4 vividly demonstrates. Such a behaviour can be easily 
explained implementing again the potential well model [53]. Introduction of the one 
Neumann disc breaks the wave guide symmetry with respect to the plane z — 1/2. This 
is reflected in the corresponding asymmetry of the wave function which inside the disc, 
for the fixed radius r < a, reaches the maximum at z = and continuously decreases 
to zero at z — 1. On the contrary, second identical Neumann window restores this 
symmetry with the function / reflecting it in figure 4. 

It was already mentioned in the Introduction that bound states are characterized 
by the two quantum numbers: principal n and azimuthal m. To exemplify this, we plot 
in figure 5 functions /' = exp(—im6) f for the fixed radius a = 3 and several quantum 
numbers (n, m) for both configurations of the boundary conditions. In addition to the 
discussed above features, one can notice that the state with the principal number n has 
ra + l maxima on r axis at the fixed z. For the nonzero angular momentum the wave 
function is zero at r = while for m = it achieves here the largest maximum. All 
these properties are in accordance with the general rules of quantum mechanics for the 
radially symmetric systems [53]. Note also that no any radial maxima and minima are 
observable in the semi infinite DD region: all of them are located in the NN and/or DN 
part, and the function / exponentially decreases to zero as the radius swipes the r axis 
to the infinity. 

4-2. Two different Neumann discs: anticrossings and their evolution 

Next, consider the case of the two concentric Neumann discs with different radii a and 
b on the opposite walls of the otherwise uniform Dirichlet channel. Only the region 
b/a < 1 is under investigation since, as stated above, we assume that b < a, b > 0. 
Here, one needs to consider three spatial regions: first, for r < b, when the transverse 
motion is confined by the two Neumann plates; second, for b < r < a when the motion 
in the z direction is governed by the Neumann requirement on the one side and the 
Dirichlet demand on the opposite walls; and, the last, the region of r > a with the pure 
Dirichlet transverse boundary conditions. Total solution of (7) for the first region is 
described by equation (20), for the region of r > a it is written as equation (196), and 
for the ring b < r < a one has: 




j=0 




(29) 
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where the Bessel functions J\ m \ (^\J E — (j + l/2) 2 r^ and Y\ m \ i^\J E — (j + l/2) 2 r^ 
for _E < (j + 1/2) 2 transform into the modified Bessel I\ m \ (^\J (j + 1/2) 2 — Er^j 

and MacDonald K\ m \ (k\J (j + 1/2) 2 — Er^j functions, respectively, and unknown 

coefficients Cj m ' and are to be found from the matching of the function / and 
its radial derivative at the circles r = b and r = a. This matching leads again to the 
eigenvalue equation of the form of (26), where, however, in the present case the infinite 
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Figure 5. Functions /' = e~ lmS f in terms of z and r for the Neumann radius a = 3 
and several quantum numbers n and m. Two capital characters denote a corresponding 
type of the boundary conditions on the opposite wave guide walls inside the circle r < a. 
Note different /'-axis scales for each of the panels. 



matrix Q takes the block form: 



Q 



Q( 3 ) QW 



(30) 



Here, the elements of the infinite matrices QW, i = 1,2,3,4, can be derived in a 
way similar to the previous subsection; in particular, in addition to the matrices 
and p( 2 \ they also contain the matrix P^ 3 ) describing the coupling between the pure 
Neumann transverse motion and its Neumann-Dirichlet counterpart: 

' iW+7'2 j'+l/2 ■ / q 

I l ) ir ,'' + 1/2 > J U " 



p(3) 



xf D (z)x? N i 



z)dz 



(31) 



7T i' + l/2' 

Energies E n0 are plotted in panels (a)-(d) of figure 6 as a function of the radii ratio for 
several fixed outer radius a. Similar pictures are obtained for the nonzero azimuthal 
numbers m too. As expected, for the vanishing inner radius, & — ?- 0, the energies 
approach those of the one Neumann window case [39] while for the b — >■ a one recovers 
the two equal discs geometry discussed in the previous subsection. Note quite large 
slope dE/d (b/a) of the lowest curves depicting the case of the large outer radii a > 10 
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b/a 




Figure 6. (a)-(d) Energies E n o as functions of the ratio b/a for several values of the 
outer radius a where the panel (a) is for a = 5, (b) - for a = 10, (c) - for a = 15, 
and (d) - for a = 20. Solid lines denote the ground state (n = 0), dotted lines depict 
the energies of the first excited level (n = 1), dashed lines are for n = 2, and dash- 
dotted lines - for n — 3. Note different energy scales for each of the figures. Arrows in 
panel (b) show energies for which the corresponding functions / are plotted in figure 
7. Numbers near the arrows denote corresponding inner radius b. Panels (e) and (f) 
depict the mean radius r„o as a function of the inner radius b for a = 10 [panel (e)] 
and a — 20 [panel (f)]. The same convention as in panels (a)-(d) is used. Different 7 
and b scales are used in each of the panels. 
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Figure 7. Wave function /^o in terms of r and z for the outer Neumann radius a = 10 
and several inner radii b. Corresponding energies £2,0 are shown by the arrows in panel 
(b) of figure 6. 



[panels (b)-(d)]. The point is that here at both limiting parameters of the inner radius 
the ground-state energy lies very close to the corresponding fundamental threshold: 
Eq d = 1/2 for 6 = and Eq n = for b = a. Accordingly, a change of the inner radius 
forces the particle to undergo transitions from the one situation to the other one and, 
since this energy interval is large compared to the smaller a, one observes rather steep 
energy descent on b/a axis. 

Another remarkable feature of the energy spectrum is emergence and evolution 
of the avoided crossings between the levels with the same azimuthal number m and 
different principal numbers n. A hint of the anticrossing formation is already seen for 
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the moderate outer radius a = 5. As this radius increases, the interval between the 
levels at b = diminishes, as it was discussed in the previous subsection and reference 
[39]. For the moderate inner radius b its change almost does not affect the energies. 
This is exemplified in panel (b) of figure 6 where the arrows show the energies at which 
the functions / are plotted in figure 7 for the second excited level, n = 2. It is seen that 
for the nth level the energy stays almost constant with b growing until the expanding 
NN disc absorbs enough of the first maximum to alter the energy E. Further growth 
of b and its subsequent approach to the first radial node drops the energy of the level 
to its lower counterpart and significantly distorts the NN maximum. For example, for 
the case of b = 0.8 in panel (b) of figure 6 the corresponding part of figure 7 shows the 
conspicuous alteration of the first extremum with the first node located at r ~ 2. As 
soon as the expanding inner NN part of the wave guide approaches the next extremum of 
the function /, the energy E nm becomes strongly 6-dependent. The rapid decrease of the 
energy in this narrow inner radius range is accompanied again by the strong modification 
of the wave function /. For example, it is seen from the comparative analysis of these 
two figures that for b ~ 1.7 the radial node of the function / is close to the boundary 
between the NN inner disk and neighbouring ND annulus of the width a — b. As a result, 
we observe an almost flat plateau on the E-b dependence. As the second extremum of 
the function / is approached by the expanding NN disc, it gets strongly distorted what is 
clearly demonstrated by the corresponding panel for b = 2.6 of figure 7. This distortion 
of the wave function is accompanied by the simultaneous decrease of the corresponding 
energy E nm until it reaches the plateau of the neighbouring lower lying level with smaller 
n pushing it downwards. This means that the function extremum fully penetrated into 
the NN disc; accordingly, the new flat part of the energy spectrum is formed when the 
boundary between the NN and ND geometries approximately coincides with the node 
of the function [see arrow for b = 3.5 in panel (b) of figure 6 and the corresponding plot 
in figure 7]. The gap between the nth and (n — l)th levels at the point of their closest 
encounter depends on the outer radius decreasing as a increases. The sharpness of the 
transition between the plateaus increases for the larger outer radius too. Similar decrease 
of the anticrossing gap was predicted to exist for the other structures, for example, for the 
N (5-potentials with the distance between them growing [54]. Such alternating plateau 
structure with the rapid change between them in the energy spectrum continues until 
all extrema of the function / are absorbed by the expanding NN disc after which the 
energy monotonically decreases to its value at b = a [see arrow for b = 5 in panel (b) 
of figure 6 and corresponding function plot in figure 7]. Accordingly, each nth level has 
(n+ 1) plateaus. The same structure as the one presented in figures 6 and 7 is observed 
for the nonzero azimuthal numbers m too where, again, a flat plateau on the E — b 
dependence corresponds to the NN-ND boundary being located in the vicinity of the 
node of the function / when the change of the inner radius has a small influence on the 
eigenvalues of the corresponding Schrodinger equation. 
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To understand these phenomena better, it is instructive to investigate the mean 
radius r according to 

Trim nm \f\ ^ nrn) ■ (32) 

Its dependence on the inner radius b is shown in panels (e) and (f) of figure 6 for 
a = 10 and a = 20, respectively. It is seen that each plateau in the energy spectrum is 
accompanied by the same flat part of the r dependence on the inner radius. Contrary 
to the energy spectrum, the transition between the plateaus is not smooth; namely, 
the increasing b at the edge of one plateau initially pushes the wave function closer 
to the origin, as panels for b = 0.8 and b = 2.6 demonstrate. This is reflected in the 
drop of the mean radius. Subsequent growth of the inner radius brings it closer to the 
next node of the wave function what leads to the formation of the new r plateau being 
located above its predecessor since the function / is distributed wider for the larger 
inner radius. The sharp drop of the mean radius after the last flat part is explained 
again by the accumulation of the wave function near the origin (see panel for b = 5 in 
figure 7). After this pronounced minimum the magnitude of f smoothly approaches its 
value of the two equal disc geometry which, quite naturally, coincides with the mean 
radius of the system with one Neumann window: 

r| b=0 =r| b=a . (33) 

Similar to the energy, the transition between the r plateaus gets sharper for the larger 
outer radii, as the comparison of panels (e) and (f) of figure 6 demonstrates. 

Anticrossings very similar to the ones discussed above were calculated for the 
magnetic quantum ring [55] when the increasing inner radius of the field-free annulus in 
the otherwise uniform magnetic field for the large fixed values of its outer counterpart 
forms the avoided crossings of the energies of the adjacent states with the same azimuthal 
quantum number m. Passage through the energy anticrossing in this case is accompanied 
by the abrupt change of the corresponding whirling persistent current of the same 
quantum mechanical state what draws the very clear parallel to our geometry with 
the rapid variations of the mean radius r. This similarity becomes almost complete if 
one recalls that the azimuthal current in the uniform magnetic field is determined by the 
expression from the right-hand side of equation (32) [55]. Energy anticrossings present 
an ubiquitous feature of the energy spectrum of the quantum systems with finite height 
potentials [56-72] and are indispensable, for example, in the description of the quantum 
Hall effect [67]. 

Finally, we provide the proof of Proposition 3.4. For the large radius, the z 
dependence of the total function f(r,9,z) is defined by xf D ( z )i equation (21a). In 
turn, general solution of the radial equation (9) for E — 1 and j = is given as 

n( \ - I C! + C 2 lnr, m = 
[r) ~ \ C ir ~\ m \+C 2 r\ m \, m^O [6 ] 

with constants C\ and C 2 - From the demand of the vanishing of the function / at infinity 
it immediately follows that C 2 = 0. In this way, the first term of the right-hand side 
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of equation (18) is obtained. Remaining two terms there can be derived from (19b) for 
j — I and j > 2, respectively, and asymptotic properties of the Bessel functions K\ m \(x) 
for the large arguments x [48,49]. We point out that the value of C\ can be fixed by, 
say, the normalization condition, equation (27). This condition requires relatively fast 

decay of the wave function at infinity in order for the integral ff® (r, 6, z) rdr to 
be convergent. As is seen from equation (18), this can be achieved for the azimuthal 
numbers \m\ > 2. For smaller \m\ the radial integral at E = 1 diverges and, accordingly, 
the corresponding state is a resonance without square-integrable function ff°\ The 
particular value of C\ in this case can be determined say, from the ^-normalization. 
We also note that the bound states and resonances at the critical radius, similar to the 
true bound states with E < 1, are degenerate with respect to the sign of the azimuthal 
number m: the two states with the opposite signs of m possess the same energy with 
their functions ff® being complex conjugate ones. 



5. Conclusions 



Our comparative analysis has confirmed that the two-Neumann-disc geometry of 
the otherwise uniform straight 3D Dirichlet wave guide, similar to its one-window 
counterpart, binds the Schrodinger particle with its energy lying below the essential 
Dirichlet spectrum. In the former case the binding is stronger than in the latter 
configuration since the two discs present a larger perturbation to the particle motion 
in the channel what results in the smaller critical values of the discs radius at which a 
new bound state emerges from the continuum. Moreover, the change of the second disc 
diameter opens up an additional channel of manipulating of the spectral and localization 
properties of the structure and leads to new phenomena such as anticrossings in the 
energy-inner radius dependence. We have also drawn parallels to the 2D geometry 
and have revealed a drastic difference between the two cases; namely, if, for the quasi- 
one-dimensional situation, the state at the critical length with its energy equal to the 
Dirichlet fundamental propagation threshold is always a resonance with non vanishing at 
infinity wave function, in three dimensions this is true for the states with zero azimuthal 
quantum number, m — 0. We have classified the state with \m\ — 1 as a resonance 
too since its function, even though vanishing at infinity, can not guarantee square 
integrability while for the larger \m\ one has true bound levels with square integrable 
function. 

In our treatment we assumed that the transverse potential is flat. Applying the 
external electric field £ perpendicular to the interfaces will turn it into the linear 
tilted one. Influence of the transverse gate voltage U = S/d on the electronic and 
optical properties of the pure Dirichlet quantum well has been the subject of intensive 
theoretical [73-78] and experimental [77] research. It was shown that the field £ 
significantly changes charged particle behaviour in these low-dimensional nanostructures 
leading, for example, to the quantum-confined Stark effect [77]. One can expect that 
for our geometry, similar to the curved 2D Dirichlet wave guides [79], the perpendicular 
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electric field will have a noticeable influence on the bound states too. Such calculations 
are now in progress. 
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